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Abstract 



We describe a general operational method that can be used in the analysis 
of fractional initial and boundary value problems with additional analytic 
conditions. As an example, we derive analytic solutions of some fractional 
generalisation of differential equations of mathematical physics. Fraction- 
ality is obtained by substituting the ordinary integer-order derivative with 
\ Jf~) , the Caputo fractional derivative. Furthermore, operational relations be- 

1/^ ' tween ordinary and fractional difTerentiation are shown and discussed in 

("^ ' detail. Finally, a last example concerns the application of the method to 

the study of a fractional Poisson process. 

^^ , Keywords: Fractional differential equations; Operational methods; Frac- 
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1 Introduction 

/\^ . Recently, great interest has been devoted to the apphcation of fractional calcu- 

^ ' lus modelling to different fields of science e.g. rheology, biology, geomorphology 

- ■ ■ (see for a recent review [jj, |2| ) . Indeed, it is well known that the introduction 

of fractional derivatives furnishes the system with a memory mechanism that is 
of great importance, in particular for some diffusive processes [3]- It is there- 
fore important to understand the advantages of using fractional derivatives in 
classical equations. Currently, several different methods are used to solve frac- 
tional differential equations, from the classical Laplace transform for linear frac- 
tional equations [3], to Adomian decomposition [4] and recently the Homotopy 
Perturbation Method P which gives approximate solutions also for nonlinear 
fractional differential equations. In this note we provide analytic solutions of a 
general class of linear fractional differential equations by operational methods. 



In some previous papers, [y, |7| the authors investigated the differential isomor- 
phism between Laguerre and ordinary derivatives. This approach permits us to 
study generahsed evolution problems (see [8|). In our framework we generalise 
some of these results to fractional differential equations. Besides that, we find 
a way to analytically solve boundary and initial value problems (BVP and IVP, 
respectively) of some interesting equations such as the space and time fractional 
diffusion. Moreover, we present some applications considering widely used ini- 
tial conditions to highlight the reliability of our method. The organization of 
this paper is the following: in Section[2]we discuss the relation between ordinary 
derivative and Caputo fractional derivative and discuss a stochastic interpreta- 
tion of fractional differentiation from an operational point of view; in Section [3] 
we state the main result regarding the solution of linear fractional differential 
equations by operational method and give some interesting examples. 

2 Relations between ordinary derivative and Ca- 
puto fractional derivative 

In [6|, the authors derived formal operational solutions of a class of partial 
differential equations involving the so-called Laguerre derivative operator 

in terms of Tricomi functions. In particular, for the initial value problem 

(2.2) 



Vtfix,t)^0,fix,t), 
f{x,0)^g{x), 



where Ox is a linear operator with constant coefficients on x and g(x) is an 
analytic function, the operational solution can be written as 

fix,t)=Co{tOx)9ix), (2.3) 

where 

is the zeroth order Tricomi function. Furthermore they pointed out the corre- 
spondence 

A -^ ^, Co{y) -^ ey, (2.5) 

and also that if h{x,t) is the classical solution, then h{x,A^^) is the solution 
to the same problem but based on Laguerre derivative, where 



Jo 



(2.6) 



In this paper, we exploit the same framework in order to arrive at solutions to 
fractional differential equations. The Caputo fractional derivative is defined as 



D'^fix) 



1 



£^fiy) 



r(m - v) Jq {x - yy-"^+ 



jdy, z^e (0,+oo), m= M, (2.7) 



and reduces to the classical derivative when u = m. It is well-known that 
the eigenfunction associated with the Caputo derivative is the Mittag-Leffler 
function 



E,{x)=Y. 



-^ T{vr + iy 



V g (0, -I-cxd), X G 



(2. 



This can be ascertained by simply writing for example the IVP problem 

/(O) = 1, 
l/'(0)=0, 



(2.9) 



and by applying the Laplace transform to both members of (j2.9p , thus obtaining 

z''C{f}iz)~z''-'^aC{fKz). (2.10) 

Now, by considering that C{f}{z) = z''^^/(z'' — a) is the Laplace transform of 
a Mittag-Leffler function we readily arrive at 



f{t)=E,{an, t>0,i^e{0,2]. 



(2.11) 



In order to complete the correspondence between classical and Caputo deriva- 
tive, first we notice that the natural choice of the Riemann-Liouville fractional 
integral (here indicated as D^'^) is not satisfactory because 



-aor"^ 






r=0 



-^ r\r(iyr + l) 

r=0 ^ ' 



(z/,l;-ar), 



(2.12) 



where 



0(7, C; a;) = X!t 



r!r(7r + C) 



(2.13) 



is a Wright function (see J9[) and not a Mittag-Leffler function. Instead, when 
V G (0, 1], we can operate the right transformation by randomising the time t 
as follows: 



or equivalently 



/•OO 



log 



E^e" 



-1/q 



(2.14) 



(2.15) 



where S is a Wright-distributed random variable, that is with probabihty density 
function 

oo (—fY 

/5(0-E irn ^nv C>0, ^e(0,l]. (2.16) 

^r!r(l-i/(r + l)) 

We can therefore write that 



Ese-"-*" = / e-°«*VH(6dC (2-17) 

^^ r!r(l -vr-u)j, ^ ^^ 



r=0 



jfi"^ ^^ r(i 



(-l/(ar))^ 



at"" ^—1 r(l -vr -v) 
where 

^-c(j/)-Ef(:^' jyeR, (2.18) 

is the generalised Mittag-Leffler function (see 9]), and where, in the last step 

we exploited formula (5.1) of [lO|. 

Alternatively, one can define the integral operator 

©r/W = r(M)i?r/(i), ^e(o,i), (2.19) 

which is a simple modification of the Riemann-Liouville fractional integral, and 
then write that 






Beside (I2.15p . we also have 

t — > Bt". (2.21) 

In practice in our analysis the operator Bj"'' replaces the integral operator A^ 
in Dattoli et al. [6]. In the ordinary case we have that Aj^^l — J^ 1 ds — t while 
m our case 



V""! = /o Ids,. = i'^r(M)/r(l + v), where 



dsr(M) 

Note that ds,. can be linked to the interpretation of fractional integration fur- 



nished by Tarasov [11| as a measure for fractal media. 



3 Analytic solutions of a class of fractional dif- 
ferential equations by operational methods 

We here seek to solve analytically a general class of linear partial fractional 
differential equations by operational methods, using the results recalled in the 
previous section. 

Theorem 3.1. Consider the following initial value problem (IVP) 

iDU{x,t)^e^f{x,t), ^e(o,i], 

\f{x,0)^g{x), 

in the half plane t > 0, with an analytic initial condition g{x) and where Ox is a 
generic linear integro- differential operator with constant coefficients acting on x, 
and which satisfies the semi-group property, i.e. Q^Qy — Qx+y The operational 
solution of equation (|3.ip is given by: 

^-^ f (rv + f 

Proof. Using spectral properties of Caputo fractional derivative, we immediately 
have that 

D''J{x,t) = D'iE,{eQx)g{x) = QxE,{eQx)9{x) = Qxf{x,t). 

D 

Corollary 3.1. Consider the following BVP: 

fD^j{x,t)^ej{x,t), 

\fiO,t)^g{t), ('-'^ 

in the half plane x > 0, with an analytic boundary condition g{t). The opera- 
tional solution of equation p. 31) is given by: 

f{x,t) = E^{x'^Q,)g{t) = E TT-TTT^W- (^.4) 

Remark 3.1. Note that Theorem \3.1\ is still valid for v € (1, 2] if the constraint 
dtf{x,0) = holds. Moreover we can state a similar remark for Corollary \3.1\ 
but with the constraint dxfiO,t) = 0. 

Example 3.1 (IVP for the time- fractional diffusion equation). Consider the 
following IVP 

(D-f{x,t)=dxxf{x,t), ^G(0,1], 
\f{x,0) = x^ iii%-. 

Considering Theorem \3.1\ the analytic solution is given by 

fut) f np + i)t-x^-'^ .3g. 



Proof. Recalling that DJT'' = T{riy + 1)^-" /T{rv + l-v), we have that 



A^ E 



r(/3 + l)t'"'a;'3-2'- 



E' 



r(;3 + l)f-''t'"'a;'3-2r 



^ r(ri/ + 1)T{I3 + 1 - 2r) j ^^ r(ri/ + 1 - v)T{P + 1 - 2r) ' 
On the other hand 



(3.7) 



dxx I E 



r(/3 + l)t'"~x'^~2r 



^ r(rj/ + l)r(/3 + 1 - 2r) I ^^ r(ri/ + l)r(/? - 1 - 2r) 



^ r(^+i)t-x^-W) 



,r-=0 



r=0 



SO with a simple change of index in the sum (r' =?' + !) we retrieve (13.71) and 
the proof is complete. 

Note that the solution p.6p is the fractional generalisation of the classical heat 
polynomials studied by Rosenbloom and Widder [12 1. 

D 

Example 3.2 (IVP for the time- fractional equation of vibrating plates [13|). 
Consider the following initial value problem for the time-fractional equation of 
vibrating plates 



'D-f{x,t) = -dxxxxf{x,t), J^e(0,2], 
f{x,Q) = sin.T, 
[dtfix,0)^0, 

The analytic solution is given by 



(3.9) 



f{x,t) = sin a; y 



{-lyf- 



^r(ri/ + l) 



sin xEi^l—t" 



(3.10) 



Note that this is simply the solution by separation of variables. 



Example 3.3 (BVP for the space- fractional diffusion equation). Consider the 
following BVP 



D-f{x,t)^dtf{x,t), ^.£(0,2], x->0, 
f{0,t) = e~\ 
[dj{0,t) = 0. 

The analytic solution is given by 



(3.11) 



f{x,t)=e-'J2 



{-lyx"''' 



= e-'E,{-x''). 



(3.12) 



It is straightforward to realise that this solution can also be derived by separation 
of variables. 

Example 3.4 (Fractional Poisson process). Consider the differential equation 
governing the probability generating function C^i^Ujt), \u\ < I, of a fractional 
Poisson process llGj N^{t), t>0, with rate k > 0; 



D'^G'^iu, t) = -k(1 - u)G''{u, t), 
G(u,0) = 1. 



(3.13) 



It may be immediately realised that the solution to (j3.13p is 

G''{u,t)=E^{-K{l-u)e), \u\<l. (3.14) 

The rem,arkahle thing here is that we can treat the related fractional difference- 
differential equations governing the state probabilities p'^{t) — Pr{N'^{t) — k}, 
k > 0, with the same operational techniques. In particular it is known thatp'^it) 
satisfy the equations 



iD'^plit) = -\plit) + Xpl^^it), V e (0, 1], t > 0, 
Where dkfl 'is the Kronecker delta, i.e. 

We can rewrite equation p.lSp as 



I D-^plit) = -A(l - B)pl{t), V e (0, 1], t > 0, 
\pm = 4,0. 



(3.15) 



(3.16) 



(3.17) 



where B is the so-called backward shift operator and is such that B{Sk.o) = 
(5fe-i,o o,nd B^{Skfl) — B'^~^{B(Skfi)) = Sk-r.o- By considering that the eigen- 
function associated with the Caputo fractional derivative is the Mittag-Leffler 
function, we obtain that the operational solution to p.l7p is 

P^(^) = E FT— 7t(1 - ^)''5,.,o (3.18) 

r=0 ^ ' h=Q ^ ' 

Z^J{^r + l){k)^ > 

Note that (j3.18p coincides with formula (1.4) of Beghin and Orsingher flJl] . For 
the specific case v = 1, that is for the classical homogeneous Poisson process, 
we can retrieve the Poisson distribution either by specialising p.lSp or by oper- 
ational methods as follows. 

p,(t)=f: tMi_a^,„ (3.19) 

r=0 

(-At) 



Ei^E I (-i)'*-." 



r! ' — ' \h 

r=0 h.=0 

^ k\r\ ^ ' 

k\ 
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